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We explore exchange coupUng of a pair of spins in a double dot and in an optical lattice. Our 
algorithm uses the frequency of exchanges in a bosonic path integral, evaluated with Monte Carlo. 
This algorithm is simple enough to be a "black box" calculator, yet gives insights into the role of 
correlation through two-particle probability densities, visualization of instantons, and pair corre- 
lation functions. We map the problem to Hubbard model and see that exchange and correlation 
renormalize the effective parameters, dramatically lowering Ur at larger separations. 
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Lattice models are popular in solid state physics, often 
serving as simple models for atomic orbitals, especially 
in the theory of magnetism [1]. Quantum dot (QD) ar- 
rays and optical lattices (OL) have emerged as new re- 
alizations of lattices. These artificial lattices are candi- 
dates for quantum computers, where spins on exchange- 
coupled QDs act as qubits for universal quantum com- 
putation [2, S] • A basic concept of lattice models is the 
intersite exchange, in which virtual hopping leads to an 
effective spin coupling of neighboring sites. A two-site 
model is one of the simplest quantum problems, yet the 
quantitative mapping from a three-dimensional model of 
a recent double QD or OL experiment to an effective 
two-site Hamiltonian has many subtleties requiring care- 
ful treatment of exchange and correlation [H, 0, S] • 

In this Letter we develop a simple path integral Monte 
Carlo (PIMC) approach for extracting accurate singlet- 
triplet splitting from a two- or three-dimensional con- 
tinuous model. This two-particle problem has been pre- 
viously solved with direct diagonalization (DD) methods 
with a careful choice of basis functions [4] and is amenable 
to variational or diffusion quantum Monte Carlo (QMC) 
[sl. However, the simple and elegant PIMC approach is 
a more direct solution without variational bias or basis- 
set issues and offers theoretical insights into this impor- 
tant problem. We first show that the splitting energy, J, 
is easily extracted from the average permutation of the 
two-particle path integral, even when J <C ksT. This 
PIMC algorithm can be a black-box calculator, provid- 
ing accurate numerical estimates of J for technologically 
relevant models of quantum dots or optical lattices, with 
arbitrary interactions and confinement potentials. The 
path integral also allows us to ask questions about the 
correlated behaviour of the electrons. For example, do 
the electrons exchange across the barrier simultaneously, 
or do they briefly double occupy the dot? Or, how does 
the motion of one electron over the barrier correlate with 



the location of the other electron? We answer these ques- 
tions by viewing representative trajectories (instantons) 
in the path integral, and, more quantitatively, by calcu- 
lating pair correlation functions. Next, we include the 
effects of a magnetic field on electron exchange in double 
QDs using a Berry's phase for the magnetic flux enclosed 
by the electron paths [6]. Finally we use the method to 
model recent experiments of exchange coupled atoms in 
an optical trap, demonstrating its broad utility 

The mapping from a continuous model with interacting 
particles to a lattice model introduces subtle complica- 
tions. For a non-interacting system it is reasonable to 
reduce the Hilbert space to include just one orbital per 
site, such as a Wannier function centered on each poten- 
tial minimum. The non-interacting many-body ground 
state is a product state of these single particle orbitals, 
and low-lying excited states are also spanned by this ba- 
sis, so an effective lattice model is an excellent approx- 
imation. Interactions are typically added to this lattice 
model, often as on-site energies, /7, or intersite terms, V. 
For small t, this gives the well-known J = -4t^/{U-V). 

There can be a serious flaw when considering inter- 
actions in this order, by first mapping to a lattice then 
adding interactions. When interactions are added to the 
continuum Hamiltonian, correlation enters as virtual ex- 
citations to higher energy orbitals. At first this seems in- 
significant, since there may be still a one-to-one mapping 
to an effective lattice model. But, when choosing effective 
lattice parameters, one must remember that many-body 
states in the continuum model have quantum fluctuations 
that are simply not present in the lattice model. 

As a specific example, consider two electrons in two 
coupled QDs. This system is often represented at a two- 
site Hubbard model, where the sites represent the Is 
ground states of the quantum dots. Correlation terms 
involve virtual excitation of the electron to the 2px. and 
2py states of the dots. These quantum fluctuations give 
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rise to van der Waals attraction, in addition to the usual 
mean-field Coulomb repulsion. Van der Waals attrac- 
tion and other correlations renormalize the interaction 
parameters to new values, Ur and Vr. 

When we consider hopping between sites, more com- 
plications emerge. The hopping barrier has contributions 
from both the external potential and electron-electron 
interactions. While the mean-field Hartree contribution 
can simply be added to the effective potential, the fiuc- 
tuating part is not so trivial. In the transition state, an 
electron passes over a barrier whose height has quantum 
fiuct nations. Thus we expect interactions to renormal- 
ize the hopping constant, t^. At the Hartree- Fock level, 
Hund-Mulliken theory already predicts a renormalized tr 
and Ur due to long-range exchange. However, the neglect 
of correlation in Hund-Mulliken theory can lead to catas- 
trophic failure at intermediate dot separations [4]. Our 
PIMC approach includes all correlation and illuminates 
its role in barrier hopping with the concept of instantons. 

We start with the two-dimensional model for the GaAs 
double quantum dot studied in Ref. 
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+ V;,t(ri) + V;xt(r2), (1) 



with m* = 0.067me and e = 12.9. The external potential 
comes from two piecewise-connected parabolic potentials, 

KxtW = |mcJo{min[(x - d)^, {x + d)^] + y^}, (2) 



with minima at x 
oscillator length tq 



We report d relative to the 
The two lowest en- 



ergy two-electron states are spatially symmetric and anti- 
symmetric under exchange, with energies and re- 
spectively. The exchange coupling, J = has been 
calculated previously using DD on a basis of Fock states 
built from seven single particles states [4] . Much care was 
taken to test convergence with the number of states and 
careful evaluation of coulomb matrix elements. We note 
that the same quality of DD calculation in three dimen- 
sions would typically take more single particle states. 

QMC techniques give essentially exact answers to 
many problems without basis set convergence issues, and 
often work just as easily in two or three dimensions. 
PIMC is nice for QD problems Q because it does not re- 
quire a trial wavefunction. However, a direct calculation 
of either or S- with PIMC often has a large statisti- 
cal errors in energy (~ 1 meV in QDs). Instead, we use 
particle exchange statistics to estimate energy difference 
J to high accuracy (~ 1 jaeV) in PIMC. 

To calculate J, we split the partition function into two 
parts that are either spatially symmetric and antisym- 
metric under exchange, Z = Z+ + Z_. These two terms 
can be expressed as a symmetrized and antisymmetrized 
imaginary-time path integral, [1, [13] 
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FIG. 1: (Color online) PIMC results, (a) Exchange couplings 
J for huj — A meV (•) and huj — Q meV (o) double QDs with a 
piecewise parabolic potential (inset). Dashed lines are direct 
diagonalization results from Ref. 4. (b) The double dot occu- 
pation probability x. Using J and x we fit (c) tr and (d) Ur 
parameters for an effective two site Hubbard-Model. Dashed 
lines in (c) show the bare hopping t for a single electron in the 
double QD, while the dashed line (d) is for V — /ed and U 
taken from a PIMC calculation on a single QD. 



This is a sum over all two-particle worldlines R{t) with 
the boundary condition R{Ph) = Pi?(0), where the oper- 
ator V swaps the particles. The symbol (±1)^ takes on 
the values (±1)^ = 1, (+1)^ = 1, and (-1)^ = -1. At 
low temperature, only the two lowest states contribute 
to the partition functions, so Z± = e~^^^ . Then, 



-/3J 



((-1)^)+. 



P=I,V 



(4) 

or J = —kBT\n{{—l)^)-^. Thus the exchange coupling 
can be calculated by sampling a symmetric (bosonic) 
path integral and taking the average of (—1)^, which 
is +1 for identity paths and —1 for exchanging paths. 

We ran PIMC simulations [lO|with our open-source pi 
code for the dots studied in Ref. [j, with the results shown 
in Fig. [T] (a). While each geometry was sampled for four 
hours on eight cores, the algorithm can be ran for just one 
minute to get quick answers with larger error bars. To aid 
other researchers, we have made the simulation available 
as a tool on nanoHUB [ill]. Coulomb interactions are 
included with a crude pair action that correctly handles 
the cusp condition. We observed convergence of the path 
integral results with 6400 discrete slices, but a higher- 
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FIG. 2: (Color online) Paths and pair densities for a double 
dot. (a) Simplified instanton with double occupation of the 
right dot. (b) Pair densities p(xi,X2) with the lowest density 
contour line that conects both potential minima (+ markers) 
at {±d,^d). (c) Simplified instanton with simultaneous ex- 
change, (d) Actual path showing double occupation, sampled 
from hco = 4 meV, d = 1.5ro QDs. (e) Actual path showing 
simultaneous exchange, sampled from hou = 6 meV, d = 2ro 
QDs. Trajectories (d) and (e) are also plotted in (b). 



quality pair action \lO\ could require fewer slices. We 
see near perfect agreement with DD, and speculate that 
small deviations are due to the finite basis in the DD 
calculation, which can be a problem as d increases 

To learn more from our simulation we collect the two- 
particle density, ^2), which is the probability to find 
one electron at xi and the other at X2, integrated over 
all values of yi and ?/2, shown in Fig.[2](b). We calculate 
double occupation, xd, which we define as the probabil- 
ity for the electrons to lie on the same side of the x = 
plane. From J and we are able to deduce renormal- 
ized values for tr and Ur — Vr in an effective two site 
model. Fig. [T] (c) and (d). Interactions renormalize tr 
to smaller values. This is consistent with Hund-Mulliken 
theory or a larger renormalized mass. The larger J comes 
from the dramatic decrease — Vr at larger dot sep- 

arations, as correlation enables more virtual hoping. 

There are two minima, (x\^X2) = T^^)^ in the 
external potential Vextij^i) + Ka;t(r2), marked in 
Fig. [2](b). For non-zero J, some paths must go between 
these minima. In a semiclassical model, the paths fluc- 
tuate around the potential minima, with rapid cross- 
ings called instantons, in which particles exchange be- 
tween the dots. An instanton can involve brief double- 
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FIG. 3: (Color online) Crossing density, p2(x,x), equivalent 
to the diagonal of the pair densities in Fig. [2jb). 



occupation of a dot, illustrated in Fig. [2ja), or simul- 
taneous exchange, as in Fig. [2](c). In Fig. [2]^d) and (e) 
we show actual paths sampled from our simulation that 
resemble the idealized instantons. In Fig. [2l^c), one in- 
stanton can be seen to move from the ((i, —d) minimum, 
briefly double-occupy the left dot, {—d^—d), then move 
to the {—d^d) minimum, while the other instanton moves 
directly between the two minima. 

Contours of X2), Fig. [21(b), illustrate a trend with 
increasing dot separation. For small d the highest prob- 
ability is directly between the minima (simultaneous ex- 
change), but at larger d the highest probability has two 
pathways (brief double occupation). Fig. [3] show the 
probability density for crossing, p(x, x). Crossing is most 
likely in the middle {x = 0) when the dots are close to- 
gether. When the dots are farther apart, the crossing 
probability has a double peak near x = id, while the 
probability for crossing in the middle is approximately 
half the peak value. The double peaks are slightly larger 
for the wider huj = 4 meV dot, consistent with more dou- 
ble occupation, but the difference is not very pronounced. 

To underscore the presence of electronic correlation 
during tunneling, we plot the correlation hole form in 
Fig.H with PIMC results next to DD results While 
some quantitative differences are apparent, consistent 
with the finite basis size in DD, the overall agreement 
is quite good. The message is clear: in the instanton, as 
one electron moves between the dots, the other electron 
moves away, enhancing the instanton and increasing J. 

For charged particles, magnetic fields can be used to 
tune the exchange coupling and even change its sign 
In the path integral, a magnetic field is easily imple- 
mented as a Berry's phase q^B, where q is the electron 
charge and <I>b is the total magnetic flux enclosed by the 
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FIG. 4: (Color online) Conditional density of one electron 
when the other electron is in the white box, showing the 
correlation hole during an instanton. Panels (a)-(c) are the 
hoo = 4 meV, d = 1.5ro QDs and (d)-(f) are the hoo = 6 meV, 
d — l.Oro QDs. Numerical factors are the likelyhood of the 
first electron being in the white box. PIMC results are shown 
on the left of each image, with DD results [4] on the right. 




FIG. 5: (Color online) Magnetic field dependence included 
with a Berry's phase for several double QDs. 



worldlines of the two electrons. The exchange splitting 
is then J{B) = -kBT\n{{e'^^^{-l)^)^/{e'^^^)^). The 
quantities are averaged from the bosonic path integral 
with no field, so data for different B-fields may be col- 
lected simultaneously. For very large magnetic fields the 
expectation value in the denominator is small and Monte 
Carlo sampling errors are catastrophic. In practice, we 
find that fields up to 4 T in strength are practical for the 
geometries we study, yielding the results in Fig. O 

As a final example, consider the exchange of two ^^Rb 
atoms in a double- well optical trap 0- This system re- 
sembles the double QD, only with much heavier parti- 
cles, a much lower temperature, short range interactions. 
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FIG. 6: (Color online) Spin-splitting of ^^Rb atoms trapped in 
a double well: x , PIMC results at 10 nK, and o, experimental 
data 7]. Insets show atomic probability densities. 



and a different confining potential. While the experi- 
ments in [7] have very little correlation, we present results 
here to show the method is practical for such systems 
and could be used for predictions of strongly interact- 
ing systems. The experiment has a double- well poten- 
tial, V{x) = Mong sin^(7rx/A) + Fghort cos^(27rx/A), with 
A = 765 nm and Mong lOEr, where Er = /i^/2MRbA^ 
f^. We model interactions as V{r) = Vq sech^ K,r with 
Vb = 50.5 /iK and hz = 0.1 nm~^ to reproduce the ^'^Rb 
scattering length. Fig [6] shows J as the barrier Fghort is 
varied, confirming agreement with experiment. 

In conclusion, we have demonstrated a PIMC algo- 
rithm for computing exchange-splitting in double QDs 
and optical lattices. The exchange splitting arises from 
instantons in the path integral, and we have collected 
statistics for these path crossings, including double occu- 
pation and the correlation hole. Correlations renormalize 
tr and Ur — Vr^ with dramatic decrease in /7r — — r at 
large separation. We find that simultaneous crossing oc- 
curs more often with closely spaced dots, while further 
separated dots are more likely to have instantons with 
double occupations. Finally, we have demonstrated the 
versatility of the algorithm with the inclusion of magnetic 
fields and applications to laser-trapped atoms. 
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